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The production of jets is studied in collisions of virtual photons, 7*p and 
7*7*, specifically for applications at HERA and LEP2. Photon flux factors 
are convoluted with matrix elements involving either direct or resolved pho- 
tons and, for the latter, with parton distributions of the photon. Special 
emphasis is put on the range of uncertainty in the modeling of the resolved 
component. The resulting model is compared with existing data. 
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1 Introduction 



The photon is a complicated object to describe. In the DIS region, i.e. when it is 
very virtual, it can be considered as devoid of any internal structure, at least to first 
approximation. In the other extreme, the total cross section for real photons is dominated 
by the resolved component of the wave function, where the photon has fluctuated into a 
qq state. The nature of this resolved component is still not well understood, especially 
not the way in which it dies out with increasing photon virtuality. This dampening is 
likely not to be a simple function of virtuality, but to depend on the physics observable 
being studied, i.e. on the combination of subprocesses singled out. 

Since our current understanding of QCD does not allow complete predictability, one 
sensible approach is to base ourselves on QCD-motivated models, where a plausible range 
of uncertainty can be explored. Hopefully comparisons with data may then help constrain 
the correct behaviour. The ultimate goal therefore clearly is to have a testable model for all 
aspects of the physics of 7*p and 7*7* collisions. As a stepping stone towards constructing 
such a framework, in this paper we explore the physics associated with the production 
of 'high-px' jets in the collision. That is, we here avoid the processes that only produce 
activity along the 7*p or 7*7* collision axis. For resolved photons this corresponds to the 
'soft' or 'low-p_i_' events of the hadronic physics analogy, for direct ones to the lowest-order 
DIS process 7*q — > q. 

The processes that we will study here instead can be exemplified by 7*7* — > qq (direct), 
7*g — > qq (single- resolved for 7*7*, direct for 7*p) and gg — > qq (double-resolved for 7*7*, 
(single-)resolved for 7*p), where the gluons come from the parton content of a resolved 
virtual photon or from the proton. Note that these are multi-scale processes, at least 
involving the virtuality Q 2 of either photon (i — 1,2) and the p\ of the hard subprocess. 
(In 7*7* physics the notation P 2 is often used instead of Q 2 , especially for the less virtual 
of the two photons; here we will use Q 2 throughout, however.) For a resolved photon, 
the relative transverse momentum k± of the initial 7* — > qq branching provides a further 
scale, at least in our framework. This plethora of scales clearly is a challenge to any 
model builder, but in principle it also offers the opportunity to explore QCD in a more 
differential fashion than is normally possible. 

At large photon virtualities, a possible strategy would be to express the cross sections 
entirely in terms of processes involving the photon directly, i.e. to include branchings 
such as 7* — > q'q' and q' — ► q'g in the Feynman graphs calculated to describe the process, 
so that e.g. the 7*p process gg — > qq is calculated as 7*g — > q'q'qq. With decreasing 
virtuality of the photon, such a fixed-order approach is increasingly deficient by its lack of 
the large logarithmic corrections generated by collinear and soft gluon emission, however. 
Furthermore, almost real photons allow long-lived 7* — > qq fluctuations, that then take on 
the properties of non-perturbative hadronic states, specifically of vector mesons such as 
the p°. It is therefore that an effective description in terms of parton distributions becomes 
necessary. Hence the resolved component of the photon, as opposed to the direct one. 

That such a subdivision is more than a technical construct is excellently illustrated 
by the x° bs plots from HERA [|lj. This variable sums up the fraction of the original 
photon light-cone momentum carried by the two highest-E'x jets. A clear two-component 
structure is visible. The peak close to x° bs = 1 can be viewed as a smeared footprint of 
the direct photon, with all the energy partaking in the hard interaction, while the broad 
spectrum at lower x° bs is consistent with the resolved photon, where much of the energy 
remains in a beam jet. The distinction between the two is not unique when higher-order 
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effects are included, but it is always possible to make a functional separation that avoids 
double-counting or gaps. 

The resolved photon can be further subdivided into low-virtuality fluctuations, which 
then are of a nonperturbative character and can be represented by a set of vector mesons, 
and high-virtuality ones that are describable by perturbative 7* — > qq branchings. The 
former is called the VM D (vector meson dominance) component and the latter the anoma- 
lous one. The parton distributions of the VMD component are unknown from first prin- 
ciples, and thus have to be based on reasonable ansatze, while the anomalous ones are 
perturbatively predictable. This separation is more ambiguous and less well tested than 
the one between direct and resolved photons. In principle, studies on the structure of 
the beam remnant, e.g. its p± distribution, should show characteristic patterns. Unfortu- 
nately, the naively expected differences are smeared by higher-order QCD corrections (es- 
pecially initial-state radiation), by the possibility of multiple parton-parton interactions, 
by hadronization effects, and so on. (Experimentally, gaps in the detector acceptance, 
e.g. for the beam pipe, is a further major worry.) Many of these areas offer interesting 
challenges in their own right; e.g. the way in which multiple interactions die out with 
virtuality, both that of the photon itself and that of the qq pair it fluctuates to. Models 
for one aspect at the time are therefore likely to be inadequate. Instead we here attempt 
a combined description of all the relevant physics topics. 

The traditional tool for handling such complex issues is the Monte Carlo approach. 
Our starting point is the model for real photons |2j and the parton distribution param- 
eterizations of real and virtual photons already present in the Pythia j|] generator. 
Several further additions and modifications have been made to model virtual photons, as 
will be described in the following. Other generators with an overlapping scope include, 
among others, HERWIG §, LDC §, LEPTO 0, PHOJET § and RAPGAP 0. The 
details of the approaches are different, however, so this gives healthy possibilities to com- 
pare and learn. Another alternative is provided by matrix-element calculations HIOH , that 
do not provide the same complete overview but can offer superior descriptions for some 
purposes. 

The plan of this paper is the following. In section 2 the model is described, with special 
emphasis on those aspects that are new compared with the corresponding description for 
real photons. Thereafter, in section 3, the dependence of simple observables on model 
parameters is illustrated. Comparisons are shown with some sets of data from ep and 
e + e~ colliders, and this is used to constrain partly the freedom in the model. Based on 
this experience, some further observables are then proposed and studied, to help shed 
light on the nature of the virtual photon. Finally, section 4 contains a summary and 
outlook. 



2 The Model 



The electromagnetic field surrounding a moving electron can be viewed as a flux of pho- 
tons. Weizsacker [|ll|] and Williams |12|] calculated the spectrum of these photons, neglect- 
ing the photon virtualities and terms involving the longitudinal polarization of photons. 
This approximation is well-known 
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to be a good approximation when the scattered 
lepton is tagged at small scattering angles. 

In the equivalent photon approximation |14| , the cross sections for the processes ep —>■ 
eX and ee — > eeX, where X is an arbitrary final state, can then be written as the 
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convolutions 



/l 
— N(u) da( 7 p -> X) (1) 

and 

dcr(ee -> eeX) = / 1 ^1^1 N{u x )N{u 2 ) da( 77 -> X) , (2) 
J J UJi u 2 

where uo is the energy of the emitted photon. In this approximation, the distribution in 
photon frequencies N{uj)duj/uj is obtained by integrating over the photon virtuality Q 2 . 
The maximum value Q^ax * s usually given by experimental conditions like anti-tagging, 
i.e. that the scattered lepton is not detected if its scattering angle is too small. 

A better approximation, and the one used in our approach, is to keep the Q 2 depen- 
dence in the photon flux f(y,Q 2 ) (with y u>/cu max , see below) and in the subprocess 
cross sections involving the virtual photon(s), 7 *p — > X and 7 * 7 * — > X, and to sum over 
the transverse and longitudinal photon polarizations. Equations ([I]) and (fj) then modify 
to 

da (ep - eX) = ^ f Uy dQ 2 f* /e (y, Q 2 ) da^v -> X) (3) 

and 

da(ee eeX ) = ffffd yi dQ 2 dy 2 dQ\ f^ c ( yi , Ql)f^ /c (y 2 , Q\) da( 7 * i7 | 2 - X) . 

(4) 

For ep events, this factorized ansatz is perfectly general, so long as azimuthal distributions 
in the final state are not studied in detail. In e + e~ events, it is not a good approximation 
when the virtualities Q 2 and Q\ of both photons become of the order of the squared 
invariant mass W 2 of the colliding photons [|15|]. In this region the cross section have terms 
that depend on the relative azimuthal angle of the scattered leptons, and the transverse 
and longitudinal polarizations are non-trivially mixed. However, these terms are of order 
QlQl/W 2 and can be neglected whenever at least one of the photons has low virtuality 
compared to W 2 . 



2.1 The Photon Flux 

When Q 2 /W 2 is small, one can derive JR], [TJ], [T5| 
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the second term, proportional to mf/Q 4 , is not leading 
log and is therefore often omitted. Clearly it is irrelevant at large Q 2 , but around the 
lower cut-off Q^in ^ significantly dampens the small-?/ rise of the first term. Overall, 
under realistic conditions, it reduces event rates by 5-10% O, [17 
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Figure 1: Schematic figure of 77 processes, illustrating the notation. 



The y variable is defined as the light-cone fraction the photon takes of the incoming 
lepton momentum. For instance, for events, Fig. [I], 

Vi = ft ' 3 = 2(1) for i = 1(2) . (7) 

Alternatively, the energy fraction the photon takes in the rest frame of the collision 
can be used, 

«»(«i + k 2 ) 

The two are simply related, 

Q 2 

Vi = Xi + — , (9) 
s 

with s — (ki + k 2 ) 2 . (Here and in the following formulae we have omitted the lepton 
and hadron mass terms when it is not of importance for the argumentation.) Since the 
Jacobian d(yi,Q 2 )/d(xi,Q 2 ) = 1, either variable would be an equally valid choice for 
covering the phase space. Small Xi values will be of less interest for us, since they lead to 
small W 2 and hence no high-p^ jet production, so y%/xi ~ 1 except in the high-Q 2 tail, 
and often the two are used interchangeably. Unless special Q 2 cuts are imposed, cross 
sections obtained with f l (x, Q ) dx rather than fl (y, Q ) dy differ only at the per mil 
level. For comparisons with experimental cuts, it is sometimes relevant to know which of 
the two is being used in an analysis. 

In the ep kinematics, the x and y definitions give that 

W 2 = xs = ys-Q 2 . (10) 

The W 2 expression for l + l~ is more complicated, especially because of the dependence on 
the relative azimuthal angle of the scattered leptons, ^12 = V 9 ! — f2- 



W 2 = x l x 2 s + ^^ -2\/l- xt-^Jl- x 2 -^Q 1 Q 2 cos <p 12 

)2n2 



Q2Q-2 

yms - Ql - Ql H — - 2a/1 - yi a/1 - 2/2<2iQ2Cosy?i2 . (11) 
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The lepton scattering angle 0j is related to Q 2 as 

$ = T^~ m * + (1 ~ Xl) ( s ~ (i-^) 2 ™' ~ 2m ") sin2 ^/ 2 ) ' ^ 

with m 2 = kf = k'1 and terms of 0(m 4 ) neglected. The kinematical limits thus are 

{QlUn « T^-"!?, (13) 

1 - a?j 

(QDmax « (l-Xi)s, (14) 

unless experimental conditions reduce the 0j ranges. 

In summary, we will allow the possibility of experimental cuts in the Xi, jji, Q 2 , 9i 
and W 2 variables. Within the allowed region, the phase space is Monte Carlo sampled 
according to Y\_AdQ 2 /Q 2 ) (^ x i/ x i) dy>t, with the remaining flux factors combined with the 
cross section factors to give the event weight used for eventual acceptance or rejection. 



2.2 Photon Processes 



The hard-scattering processes are classified according to whether one or both photons are 
resolved. For completeness we here quote some of the less familiar cross sections. 
For the direct process 7*7* — > ff, f some fermion, the cross sections are [18| 
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(15) 

,(16) 
(17) 
(18) 

(19) 
(20) 



ef the electrical charge and Nf the colour factor, 3 for a quark and 1 for a lepton. Re- 
member that s + t + u = —Qi — Q\i neglecting the fermion mass. Note that the cross 
section for a longitudinal photon vanishes as Qf in the limit Q 2 — ► 0. 

For a resolved photon, the photon virtuality scale is included in the arguments of the 
parton distribution but, in the spirit of the parton model, the virtuality of the parton 
inside the photon is not included in the matrix elements. Neither is the possibility of 
the partons being in longitudinally polarized photons (see below, however). The same 
subprocess cross sections can therefore be used for direct 7*p processes and for single- 
resolved 7*7* ones. For 7*q — > gq one obtains (Tj| 



and for 7*g — ► qq 



da T 2 1 i 2 + u 2 \ 2Q 2 s 



dt q (s + Qf) 2 £tt 



(s + Ql) 2 



(23) 



~dT = (24) 

Convolution with parton distributions gives 

da( 7 *p - X) = J J dx 2 dtff(x 2 ,fx 2 ) ^(s = x 2 W 2 ) , (25) 

da(7*y-X) = jjdx 2 diff(x 2 ^ 2 ,Q 2 2 )^(s = x 2 W 2 ), (26) 



where /x 2 is the scale of the hard-scattering subprocess. 

Finally we come to resolved processes in 7*p and doubly-resolved ones in 7*7*. There 
are six basic QCD cross sections, qq' — ► qq', qq — > q'q', qq — > gg, qg — > qg, gg — > gg 
and gg — ► qq. Since again parton virtualities are not included, these are the expressions 
familiar from pp physics |2(J and are not listed here. Again a convolution with parton 
distributions is necessary, 

da( 7 *p^x) = JJJ dx 1 dx 2 d^77*(£l^^Ql)//(£2,/i 2 )^(s = x 1 x 2 ^y 2 ) , (27) 

da( 7 V-X) =JjJ dx 1 dx 2 diff(x l ,fi 2 ,Q 2 1 )fj\x 2 ,fi 2 ,Q 2 2 )^(rs = x 1 x 2 W 2 ) .(28) 

In line with the neglect of parton masses, also the s expression for 7*7* is simpler than 
its W 2 analogue in eq. (O). When initial-state radiation is subsequently included, both 
transverse momenta and spacelike parton virtualities are generated, but in such a way that 
the relation s = X\X 2 W 2 is maintained JO]]. Differences between alternative s definitions 
being subleading, they are beyond our standard QCD accuracy. 



2.3 Parton Distributions 

One major element of model dependence enters via the choice of parton distributions 
for a resolved virtual photon. These distributions contain a hadronic component that is 
not perturbatively calculable. It is therefore necessary to parameterize the solution with 
input from experimental data, which mainly is available for (almost) real photons. In 
the following we will use the SaS distributions [BJ, which are the ones best suited for our 
formalism. Another set of distributions is provided by GRS ||22|| , while a simpler recipe 
for suppression factors relative to real photons has been proposed by DG |23|j . 
The SaS distributions for a real photon can be written as 

Pa(^ A* 2 ) = E 4 -^Pa' V ^ Ql) + E 2 < f 2 TT ^' **) • ( 29 ) 

V JV q K 

Here the sum is over a set of vector mesons V = p ,u,(p,J/ip according to a vector- 
meson-dominance ansatz for low-virtuality fluctuations of the photon, with experimentally 
determined couplings 4na cm / f v . The higher-virtuality, perturbative, fluctuations are 
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represented by an integral over the virtuality k 2 and a sum over quark species. We will 
refer to the first part as the VMD one and the second as the anomalous one. Each 
component f 1,v and f7><& obeys a unit momentum sum rule, and also obeys normal QCD 
evolution equations. The p' v (x, Qq) have to be determined by a tune to F^(x, ji 2 ) data, 
while /T'qq evolve from the boundary condition 

«(*, k 2 ; k 2 ) = ^{x 2 + (l- xf) (5 aq + 5 m ) . (30) 

The fi 2 dependence enter both via the evolution of each component and via the upper limit 
of the dk 2 integral. It is the latter dependence that generates the so-called anomalous 
term of the photon distribution evolution equations, from which the terminology has been 
taken over for the related event class. 

From the above ansatz, the extension to a virtual photon is given by the introduction 
of a dipole dampening factor for each component, 




Thus, with increasing Q 2 , the VMD components die away faster than the anomalous ones, 
and within the latter the low-/c 2 ones faster than the high-A; 2 ones. In the VMD component, 
the effective evolution range is reduced by the introduction of a Ql = Q 2 lt = QoQ c s > Qo, 
with 

Q 2 - a 2 Q ° + Q2 exp [ V ~ g°) ] (32) 
QcS /i 2 + Q 2 6XP I (/i 2 + Q 2 )(Q 2 + Q 2 )) ■ {d2} 

As a technical trick, the handling of the k 2 integral is made more tractable by replacing the 
dipole factor by a /^-independent multiplicative factor and an increased lower limit Qf nt of 
the integral, in such a way that both the momentum sum and the average evolution range 
is unchanged. Finally, correction factors are introduced to ensure that /J* (x, fi 2 , Q 2 ) — > 
for fi 2 — > Q 2 : in the region Q 2 > fi 2 a fixed-order perturbative description is more 
appropriate than the leading-log description in terms of a resolved photon. We then 
arrive at the so-called modified P mt scheme, which is the one used here. 

Since the probed real photon is purely transverse, the above ansatz does not address 
the issue of parton distributions of the longitudinal virtual photons. One could imagine 
an ansatz based on longitudinally polarized vector mesons, and branchings 7£ — > qq, but 
currently no parameterization exists along these lines. We will therefore content ourselves 
by exploring alternatives based on applying simple multiplicative factors R to the results 
obtained for a resolved transverse photon. As usual, processes involving longitudinal 
photons should vanish in the limit Q 2 — > 0. To study two extremes, the region with a 
linear rise in Q 2 is defined either by Q 2 < [i 2 or by Q 2 < m 2 , where the former represents 
the perturbative and the latter some non-perturbative scale. Also the high-Q 2 limit is 
not well constrained; we will compare two different alternatives, one with an asymptotic 
fall-off like 1/Q 2 and another which approaches a constant ratio, both with respect to the 
transverse resolved photon. (Since we put /J*(x, fj, 2 , Q 2 ) = for Q 2 > fj 2 , the R value 
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will actually not be used for large Q 2 , so the choice is not so crucial.) We therefore study 
the alternative ansatze 



fl.teC/O = H.^^, (33) 

o , ^2 2, , 4 2 /l/lfe.C 2 ) ,,,, 

w.m = H V+ w/|(»,Qr (34) 

„, n2 2, , , "<y / 7 L /<(y.Q 2 ) 

H.to.0,/.) = 1 + a K + Q2 )/ T ife , < 22 ) (35) 

with a = 1 as main contrast to the default a = 0. The y dependence compensates for 
the difference in photon flux between transverse and longitudinal photons. Ri and R 2 
have the same onset at low Q 2 but different asymptotic behaviour; for the former case the 
longitudinal part vanishes and for the latter it approaches a constant (w.r.t. the transverse 
case). For the third case the onset is governed by a non-perturbative parameter m p and 
it has the same asymptotic limit as R 2 - In a more sophisticated treatment, presumably 
also the k 2 scale of the 7* — > qq fluctuation would enter. In double-resolved 7*7* events 
one R factor is applied for each side. 

Another ambiguity is the choice of fi 2 scale in parton distributions. For a process such 
as 77 — > qq, with real photons, conventional wisdom is that /i 2 = —t is the proper scale in 
the limit t — > 0, where i-channel graphs dominate the cross sections, and fi 2 = —u in the 
limit u — > 0. The combination fi 2 = tu/s = p 2 ^ interpolates between these limits and thus 
is a traditional choice, sometimes multiplied by some constant factor. When the incoming 
(or outgoing) photons/partons are not massless, tu / s 7^ p 2 ± . A possible generalization for 
a direct virtual photon, Q\ 7^ 0, is 

2 iu tu 2 s + Q\ ( 
t + u s + Q( s 

For a corresponding t-channel graph involving the quark from a resolved photon, in princi- 
ple the same relation should hold, with the quark virtuality Q\ substituting for the photon 
Q\ one. Inside Pythia, however, all incoming partons are assumed massless in the selec- 
tion of hard-scattering kinematics. When, later on, the 7* — > qq branching is included, 
the s and the rest-frame scattering angle are left unaffected. Thus also p\ is unchanged, 
while t and u are not. Using the corrected values, one obtains the same final expression 
as in eq. (|36D , with Q\ instead of Q\. Kinematics provides the constraint Q\ > X\Q\ at 
the 7* — > qq branching. The lower bound is a rather conservative estimate of the actual 
value of Q 2 , however, and more typically one obtains Q\ ~ Q\ from the 7* — > qq branch- 
ing. (Further branchings are included in the standard QCD parton-shower description. 
By analogy with the description of other QCD processes, we here disregard these shower 
virtualities.) A scale choice like in eq. ( p6|) therefore should be a sensible one both for 
direct and resolved photons. 

When both incoming photons are virtual, the relation between /x 2 and p\ becomes more 
complicated, but in the limit that terms of order Q\Q\js 2 are neglected, they simplify to 

^ = p\ S + ( ^ + & . (37) 
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This expression does not guarantee that fi 2 > Q\ + Q\. Sometimes such an inequality 



is assumed [24], so in order to cover a broader range of scale choices, below we will be 
comparing six different alternatives. In (almost) increasing order these are 

li\ = pi, (38) 

III = pi § + XlQ l + X2Ql , (39) 
s 

l4 = P \ S + Q Y Ql , (40) 

Ha = P± + 2 ' ^ ' 

f4 = pl + Qi + Ql, (42) 

111 = 2i4. (43) 

Only the fifth alternative ensures f2* {x, fi 2 ,Q 2 ) > for arbitrarily large Q 2 ; in all other 
alternatives the resolved contribution (at fixed p±) vanish above some Q 2 scale. The last 
alternative exploits the well-known freedom of including some multiplicative factor in any 
(leading-order) scale choice. When nothing is mentioned explicitly below, the choice n\ is 
used. We should note that the expressions in the program also contain a dependence on 
final-state masses, e.g. for the production of massive quarks, but this is left out here since 
it is not a topic studied in this paper. (It will enter briefly in the following, however.) 

2.4 Other Model Aspects 

The issues discussed above are the main ones that distinguish the description of processes 
involving virtual photons from those induced by real photons or by hadrons in general. 
In common for the tree is the need to consider the buildup of more complicated partonic 
configurations from the lowest-order 'skeletons' defined above, (i) by parton showers, (ii) 
by multiple parton-parton interactions and beam remnants, where applicable, and (Hi) by 
the subsequent transformation of these partons into the observable hadrons. The latter, 



hadronization stage can be described by the standard string fragmentation framework |25| . 
followed by the decays of unstable primary hadrons, and is not further discussed here. In 
the following we comment on the shower, multiple-interaction and beam-remnant aspects. 

The parton-shower description is conveniently subdivided into initial- and final-state 
radiation. For wide-angle emissions such a classification is not unambiguous, and it is 
necessary to consider interference effects in order to describe the data [pq| . A direct 
photon is not associated with any initial-state QCD radiation. 

In the hadronic environment, initial-state radiation normally means an evolution of a 
spacelike branch of partons, from an initially vanishing virtuality up to the scale of the 
hard process. The 'backwards evolution' strategy |2]J allows this cascade to be recon- 
structed in reverse order. That is, from a parton b coming in to the hard scattering, the 
branching a — *> be that produced b is first reconstructed, thereafter the branching that 
produced a, and so on, down to the lower cutoff scale Qff 1 , in practice of the order of 
1 GeV. Inclusively, the effect of parton-shower histories is already taken into account by 
the fi 2 dependence of the parton distributions used to select the hard scattering, so what 
this procedure does is to associate an exclusive set of initial-state-radiation partons to 
each hard scattering. 

For a real photon, the VMD part is assumed to behave like a hadron, while the lower 
parton-shower cut-off Qg h has to be considered further for the anomalous component. In 
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the spirit of the ansatz for parton distributions, eq. the k 2 value of the 7 —>■ qq 

branching is distributed like dk 2 /k 2 . In principle, for a given physical process, there 
would be a further 'trigger bias' effect to this distribution: at large (small) x values 
a large (small) k 2 would be favoured since it would be associated with a small (large) 
evolution range. As a first approximation, this bias is neglected here, i.e. k 2 is picked flat 
in a logarithmic scale between Qq and fi 2 . When this k is larger than the default value 
of QoN the shower cut-off is increased to instead be given by k. Thus there is no shower 
evolution below the k scale of the 7 — > qq branching. Thereby the amount of collinear 
emission along the incoming photon direction is reduced for high-virtuality fluctuations. 

As a technical aside, we note that the default value Qo = 0.6 GeV is smaller than 
Qo 1 — 1 GeV. First, it should be made clear that k 2 and the Q 2 virtuality variable of 
the spacelike shower evolution are not defined in precisely the same way and therefore 
very well can differ by factors of order unity, so that also the lower cut-offs may differ. 
However, also disregarding this, there is no reason why the two cut-offs should agree. 
Whereas Qo is severely constrained e.g. by data, Q s Q h is just a crude estimate of where 
emitted partons in the cascade are so collinear that their effect may be safely neglected. 
That is, we do not exclude the possibility of parton emissions in the range between Q 
and Q s Q h , but also do not believe that their inclusion or not is essential for a description of 
final-state properties. The effects of an increased shower cut-off only show up for 7 — > qq 
branchings with k ^ 2 GeV, in terms of a reduced parton emission rate along the incoming 
photon direction. 

The generalization to a virtual photon is now straightforward, given the extension of 
fa to f2*, cf. eqs. ( P§D and (|5T|). For the VMD piece, the shower cut-off is selected to be 
the larger of Q int and Q s Q h . For the anomalous piece, a k 2 is picked flat in a logarithmic 
scale between Qf nt and /i 2 , and the shower cut-off is the larger of this k and Q s h . 

Final-state radiation from the scattered partons of the hard interaction follow the 
same pattern as in hadron-hadron collisions, and so need not be discussed here. Also 
time-like partons on the 'side branches' of an initial-state cascade can undergo final-state 
evolution, as part of the conventional shower description. The new aspect concerns the q 
'beam remnant' parton of an anomalous 7* — > qq branching, where the q is the initiator 
of the spacelike cascade. This parton can undergo a shower evolution from the 7* — > qq 
branching scale k down to the timelike shower cut-off m s h w 1 GeV, provided k is above 
the cut-off. For the VMD part of the photon, the beam remnant is assumed not to radiate. 

In general, a beam remnant contains the flavours 'left behind' when one parton initiates 
the spacelike shower that leads up to the hard interaction. Momentum conservation gives 
its kinematics. This includes both a longitudinal momentum fraction 1 — x, if the initiator 
takes x, and a transverse 'primordial recoil. For the VMD components of the photon 
this transverse momentum should be of a nonperturbative character and small, ~ 0.5 GeV, 
although there are indications of much larger values in data ||27[| . A Gaussian ansatz is 
used to pick a k± vector. For the anomalous part of the photon, we associate its k± with 
the k scale that already figured prominently above. 

Multiple parton-parton interactions p8| could be viewed as a further sophistication 
of the beam-remnant description. The basic idea is that hadronic states contain many 
partons and that therefore several perturbative interactions may occur in a hadron-hadron 
collision. Usually (at most) one of these give rise to visible high-p^ jets, while the rest only 
add to the underlying event activity. This additional component is vitally needed in order 
to reproduce the observed multiplicity and E± flow in high-energy pp collisions, at least if 
string hadronization is supposed to be universal, i.e. have its parameters constrained by 
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e + e~ data. The main unknown parameter in this approach is p± l min ~ 2 GeV, the scale 
below which perturbation theory is assumed not to be applied anymore. Physically, this 
scale can be viewed as related to an effective colour screening phenomenon: gluons with 
a wavelength larger than 1/p, min do not resolve the colour structure inside the hadron 
wave function and therefore decouple. Of course, this screening should set in gradually, 
so PjTmin * s on ly an effective parameter. To first approximation, the average number of 
interactions per event is the ratio of the jet cross section above p^min t° the total inelastic 
cross section. By now there exist convincing evidence in favour of multiple interactions 
also in the HERA data |29 . 



Almost by definition, direct and single-resolved processes do not lead to multiple in- 
teractions. Furthermore, within the spectrum of resolved components of a real photon, 
one would expect a close-to hadronlike behaviour for the VMD part and then a gradual 
fading-away of multiple interactions as the k 2 scale of an anomalous fluctuation is in- 
creased. This may partly be seen as a consequence of that the number of partons (above 
some fixed xq) decrease with increasing k 2 . Furthermore, the physical size of an anoma- 
lous fluctuation should scale like 1/k, so the screening argument above would lead to a 
scaling-up of p± l min roughly by a factor k/Q . This in part would be compensated by an 
expected decrease in the total cross section of an anomalous fluctuation, about like 1/k 2 
based on the size scaling argument, but in total there should still be a rather rapid fall-off 
with k. As a first approximation, we follow the route adopted in ||, namely to allow 
multiple interactions for VMD states to the same extent as for hadrons, and not at all for 
the anomalous states. 

When generalizing to virtual photons, the same arguments as used above for the 
anomalous component would suggest a fall-off of multiple interactions in the VMD com- 
ponent with increasing Q 2 scale. As a first guess, we have here chosen to scale up p^mm 
\/l + Q 2 /rriy, while the normalizing total cross section is scaled down like m v / [niy+Q 2 ) 2 
in accordance with the dipole ansatz. 



3 Results 

3.1 Basic Distributions 
3.1.1 7*p Partonic Cross Sections 

The contribution to the <7 7 * p parton cross section from the different direct components 
in equations (p])- (|2^) , the boson-gluon fusion and the QCD Compton processes, are 
shown in Fig. ^| as a function of the photon virtuality, Q 2 . The partons produced in the 
2^2 hard scattering subprocess are restricted to have transverse momenta p±_ larger 
than 5 GeV and the hadronic centre of mass energy W = v /s 7 » p is equal to 200 GeV. 

In the limit Q 2 —>■ the cross section vanishes for the processes involving a longitudinal 
photon. In the same limit the boson-gluon fusion process with a transverse photon, 
7tS — > QQ; is more important than the QCD Compton counterpart, 7^q — > gq; this 
kinematic region permits the gluon/quark from the proton to have a small fraction x of 
the proton momentum, x = (s + Q 2 )/(W 2 + Q 2 ). At large Q 2 , the opposite is true; x 
is large and the QCD Compton process, involving an initial quark from the proton, is 
more important than the boson-gluon fusion process. Both processes are suppressed by 
the l/(s + Q 2 ) 2 factor that provides the main dampening at large Q 2 . The longitudinal 
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Figure 2: Contributions to the <t 7 * p parton cross section from the subprocesses in 
equations (0])-(]2§), pf rton > 5 GeV and W rp = 200 GeV. The circles and the solid 
line is the sum of the transverse (T) and the transverse plus the longitudinal (T+L) 
direct processes respectively. 

processes, 7£q — > gq and 7£g — > qq have maximal cross sections when Q 2 ~ s, since 
s > 4p_|_ min this occurs for Q 2 around 100 GeV 2 . 

The proton structure function used to produce Fig. |2| was GRV 94 LO, which is the 
default in Pythia. A different choice of structure function would only change the details 
but not the conclusions made from these results. In the following the GRV 94 LO will be 
used as the proton parton distribution except where otherwise stated. 

When studying jet production with low photon virtualities it is likely that the photon 
is resolved. The photon fluctuates then into a quark- ant iquark pair that develops into a 
multiparton state and finally one of these partons scatters off a parton from the proton. 
The virtual photon parton distributions used for modeling the resolved components are 



the SaS ID and 2D distributions. The dipole factors in eq. ([31]) will dampen the resolved 
components with increasing photon virtualities. 

The relative importance of the direct components compared to the resolved compo- 
nents — the VMD and anomalous fluctuations — are shown in Fig. [3[ The direct con- 
tribution is the sum of the transverse QCD Compton and boson-gluon fusion processes 
studied in Fig. |]. As expected, the VMD components die away faster than the anomalous 
ones which in turn die away faster than the sum of the direct components. Here we also 
see the difference between the SaS ID and SaS 2D distributions, with SaS 2D having a 
higher cut-off for the separation between the VMD component and the anomalous com- 
ponent. This is reflected by the larger contribution from the SaS 2D VMD component 
compared to the SaS ID one, and vice versa for the anomalous component. 

Longitudinal resolved components are simulated with an extra factor of Ri(y, Q 2 , /x 2 ), 
eq. Q5BD, for the transverse resolved components. In Fig. [5] an a-parameter equal to 1 
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Figure 3: Contributions to the <t 7 * p parton cross section from events involving a 
direct, a VMD and an anomalous photon, p^ ar on > 5 GeV and v /s 7 * p = 200 GeV. 
To the left, the SaS ID and SaS 2D parton distributions was used, to the right 
SaS ID was used and a longitudinal resolved photon is modeled with a = 1 and 
y = 0.44 (see text). 

has been used together with the SaS ID photon parton distribution. The y value is set 
equal to 0.44 with HERA energies in mind. A noticeable difference is seen in the total 
jet cross section for Q 2 around 10 GeV 2 , since, in this ansatz, the maximal contribution 
for the longitudinal component is obtained at Q 2 = fi 2 , and the direct and the anomalous 
component is of the same order in this region. With a = 1 the sum of the resolved 
components can be almost a factor of two larger than with a = 0, the case with pure 
transverse resolved photons. There is nothing special with the choice a = 1 except that it 
models a longitudinal resolved component that is of the same order of magnitude as the 
transverse component at intermediate Q 2 values; this is the same behaviour as seen for 
the direct components. 

3.1.2 7*7* Partonic Cross Sections 

In Fig. f| the processes 7t,l7t,l — > °l°o equations (^5|)-(|l8D , have been generated at a 
centre of mass energy W = v /s 7 * 7 * of 100 GeV. For simplicity, one of the photons was 
kept at a fixed virtuality of 0.1 GeV 2 or 10 GeV 2 , respectively. The transverse momenta 
for the partons in the hard scattering was restricted to be larger than 5 GeV. 

Clearly, the dominant process is 7^7^ — > qq but, as the virtuality Q\ increases, the 
contribution from the longitudinal processes becomes more and more important. From 
eq. fll5D-([TED we notice that the 7^72 > 72 7r and 72 72 direct cross sections have a maximum 
at Q\ ~ s = s 7 * 7 » (for Q 2 <C s). The 7272 process is bounded to never be larger than twice 
the 7^72 or 72 It process. The 72 7t an d 7272 processes clearly go to zero as Q\ — > 0, the 
7t72 process remains finite and approaches the expression for the 72g matrix element, 
eq. ([Z4p, with appropriate factors exchanged (and Q\ Qf). The Q\ proportionality for 
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Figure 4: Contributions to the <t 7 * 7 » parton cross section from the subprocesses 
in equations (0)-(0), pf rton > 5 GeV and = 100 GeV. The Q\ of the 

first photon is varied while the second photon is kept at the fixed virtualities Q\ = 
0.1 GeV 2 or Q\ = 10 GeV 2 , respectively. 



the "fxll an d matrix elements gives the factor of 100 difference between the two 

cases in Fig. |[ As for the case with direct processes in 7*p, the 7^7^ — > qq processes 
(£i =T, L) are suppressed by the 1/A 2 factor for high Qf. 

With the possibility of having resolved photons, nine event classes is obtained. They 
are illustrated in Fig. || for the SaS ID and SaS 2D parton distributions, with the sec- 
ond photon having the fixed virtuality Q\ = 1 GeV 2 . Generally, as Q\ increases, the 
VMD components (of the first photon) drops first, then the anomalous components and 
at very high Q\ remains the direct components (direct-direct, direct-VMD and direct- 
anomalous). Particularly, with the choice Q\ = \ GeV 2 , the anomalous-anomalous events 
dominates at low Q\ for the SaS ID case, whereas for the SaS 2D case, the VMD compo- 
nents dominate. In this ansatz, the Q- int scheme, the parton distribution for the anomalous 
component /7*(x, /i 2 , Q 2 ) = when /z 2 < Q 2 . 

With the y-values taken equal to 0.5 (s 7 * 7 * ~ yi^Se+e-, i/s e + e - — 200 GeV) and 
a = 1 for Ri(y, Q 2 , /x 2 ), the longitudinally resolved photon contribution can be obtained 
as the difference between the Total T+L and the Total T in Fig. |5|. Again, we see that 
the contribution is most pronounced where p\ min ~ Q 2 . For this special case, it happens 
to be where the anomalous component (of the first photon) is comparable with the direct 
one. 

3.1.3 x 7 Distributions 

In ep collisions the variable x° bs is defined as 

^ = E ^\ + E ^ m (44) 



14 



SaS ID, Qj=l GeV 2 



SaS 2D, o|=l GeV 2 




Figure 5: Contributions to the <r 7 * 7 * parton cross section from the nine different 
event classes, pP arton > 5 GeV and Js^f = 100 GeV. The Q\ variation is for the 
first photon; the second photon is kept at fixed virtuality Q| = 1 GeV 2 . 



where E± ; i is the transverse energy and rji is the rapidity of the two jets with the highest 
transverse energy. E 1 is the energy of the photon. x° bs is then the fraction of the photon 
energy (or better, light-cone momentum) that goes into the production of the two highest 
transverse energy jets. In Fig. ^ the x° bs distribution is shown for the different photon 
components. A cone jet algorithm with cone radius R — 1 is used for jet finding. The 
partons in the hard scattering are restricted to have transverse momenta larger than 3 GeV 
and the jets to have a transverse energy larger than 6 GeV. Note that p^min < Fxnun ^ as 
been used since migration from pj_ < p 1 * rriin is likely to happen due to various effects; this 
is discussed further in section |3.2| . The x° bs distribution was generated for two different 



Q 2 intervals; 0.9-1.1 GeV 2 and 9-11 GeV 2 . In both cases, an x° bs value of 0.7-0.8 will 
separate most of the direct events from the resolved events. In the low-Q 2 case, the SaS 2D 
parton distribution gives less resolved events at large compared to SaS ID. This 

is because the VMD component, which is more important in SaS 2D, is dampened faster 
with increasing Q 2 than the anomalous one. For the high-Q 2 case, this difference is not 
there, since the two anomalous components are approaching each other at high Q 2 , and 
the VMD pieces are vanishingly small. 



In e + e collisions the variables x+ and x 1 are defined as JSC 



particles I ^ 31 V 



where p z is the momentum component along the z-axis of the e + e~ collision and E is 
the energy of the jets or particles. The sum over jets runs over all jets in the event 
but is limited to the two hardest jets for dijet studies. The different direct and resolved 
components are shown in Fig. [?| for Q 2 < 0.8 GeV 2 ; each combination of components 
will not be shown here. For these symmetric cuts the x^ distributions are identical, the 
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Figure 6: The x° bs distribution for the different photon components. p p ± T on > 
3 GeV, E?* > 6 GeV, 0.01 < y < 0.99 and v /s^ ~ 300 GeV. 

direct-resolved events for one case corresponds to the resolved-direct for the other. As 
expected, direct events are concentrated at large x^, direct -resolved (resolved-direct) are 
compatible with direct events in ep, and double-resolved are concentrated at small x^ 
- suggesting a cut at low x+ and x~ to separate double-resolved from single-resolved 
events. 

The SaS 2D distribution gives, with Q\ 2 < 0.8 GeV 2 , a large double-resolved contri- 
bution, dominated by the VMD-VMD events. Similarly, direct- VMD events dominate 
the single-resolved contribution. The effective parameter p^ l IDin sets the amount of mul- 
tiple interactions for the VMD component. The results shown in this section are sensitive 
to the choice of p^ min , here chosen to 1.4 GeV, which is the default in Pythia. 
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Figure 7: The x+ distribution for the direct and resolved photon components. 

^parton = ^ ^jet > 3 Q < y ^ < Q 99^ Q2^ < q g Ge y2 and 

y/s e + e - = 133 GeV. 
3.1.4 Scale Choice in Parton Distributions 

To show the dependence on the choice of scale /x 2 in parton distributions we will study ep 
events, with direct and resolved events treated separately. Three extreme choices of scales 
lA — v\i £*i = P± + Q 2 an d lA ~ 2p± are compared. For simplicity, we stay with the p± 
distribution of the hard scattering 2 — > 2 subprocess, i.e. avoid hadronization, multiple 
interactions, initial- and final-state bremsstrahlung and jet clustering effects. 

In Fig. H, the da ep /dlog 10 (jo_i_/(lGeV)) 2 distribution for the different scale choices are 
compared at two different photon virtualities, Q 2 ~ 1 GeV 2 and Q 2 ~ 10 GeV 2 . For the 
low-Q 2 case, their differences are most pronounced for resolved events and the results are 
in decreasing order in accordance with the scale choice; /xg, /is and Here the anomalous 
component is the dominant one, explaining the difference between the SaS ID and the 
SaS 2D distributions. For direct events, only a rather mild scale-breaking of the proton 
distributions enters. Since small-x values dominate, where distributions increase with /i 2 , 
the hq scale should give the largest result but the difference from the other two are within 
errors. 

For Q 2 w 10 GeV 2 , direct events dominate. The resolved results can be divided into 
two regions; p]_ < Q 2 and p 2 ± > Q 2 . Remember that the /is scale imply that the resolved 
component will not vanish even when p 2 ± < Q 2 as can be seen in Fig. || With the fi^ 
scale, the evolution range decreases with p 2 ± and finally suppresses resolved photons. The 
\i\ scale does not allow resolved photons when p\ < Q 2 . (The tiny tail of such events 
comes from charm production, where actually /i 2 = p\ + m 2 is used rather than only 
p\.) The results with the /is and ^ scales are, of course, in agreement when p\ w Q 2 
and in the tail of the p± distribution the p,Q scale gives the bigger cross section. At large 
Q 2 values, the difference between the two photon parton distributions are reduced. For 
direct events the /is scale is larger than the /xg one &t low p\ but smaller at high p 2 ± ; a 
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Figure 8: The dcr ep /dlog 10 (p^/(lGeV)) 2 cross section for different /ij scales with 
Q 2 w 1 GeV 2 and Q 2 w 10 GeV 2 , ^ ~ 300 GeV. 



corresponding pattern is expected to be seen in the cross section (again, small-x values 
dominate). However, differences are small. 



3.2 Comparisons with Data 

In this section the model is compared with data. We will not make a detailed analy- 
sis of experimental results but use it to point out model dependences and to constrain 
some model parameters. Where applicable, the HzTool |32[ routines will be used for the 
comparison with data. 

2^2 parton interactions normally give rise to 2-jet events. In leading-order QCD, 



18 



the jets are balanced in transverse momenta in the centre of mass frame of the 7*p/7*7* 
subsystem. Various effects, such as primordial k±, initial- and final-state bremsstrahlung, 
etc., tend to spoil this picture. This increases the da/dp± spectrum at any fixed 
since jets can be boosted up from lower p±. It is here important to remember that 
da/dp± is dropping steeply with increasing p±, since the matrix elements described in 
section |2~2] are divergent for pj_ — > 0. Even a symmetric smearing of jet p±, i.e. with shifts 
+Sp± and — Sp± equally probable, would thus increase da/dp±: the jets shifted upwards 
would migrate into a region less populated than those shifted downwards, and therefore 
proportionately have a larger impact. Additionally, the smearing is not symmetric, but 
normally shifts one jet upwards while the other may be shifted in either direction. (This 
holds for primordial kj_ and initial-state radiation, while final-state radiation tends to 
shift both jets downwards.) In the generation procedure, the divergences are avoided by 
the introduction of a lower cut, P^min > below which parton scatterings are not considered. 
In order to study jets above some p'xmim a P±min < Pimin ^ s rec L urr ed. Stable results are 
obtained when a hard scattering with pP arton = p^ a ^°° for the incoming partons do not 
produce any jets with p± above Pj^min- Typically p^min = §F±min or ^ ess * s required. 



3.2.1 Dijets at HERA 

Dijet angular distributions in photoproduction and deep inelastic scattering have been 



studied at HERA ||31|| . In the center of mass system of the two highest transverse energy 
jets (which are restricted to have a transverse energy above a certain value) the angle 
6* between the jets and the beam axis is expected to be distributed differently depend- 
ing on whether or not the photon is resolved. The leading-order direct QCD graphs, 
boson-gluon fusion and QCD Compton, have spin-1/2 quark propagators leading to an 
angular dependence proportional to 1/t oc (1 — |cos#*|) _1 . (The jets are generally not 
distinguishable so the absolute value of cos#* cannot be measured.) On the contrary, 
the resolved processes are dominated by Rutherford type scatterings involving a t-channel 
gluon: qq' — > qq', qg — > qg and gg — > gg. These thus give an angular dependence pro- 
portional to 1/t 2 oc (1 — | cos#*|)~ 2 i.e. the dijet angular distribution is expected to rise 
steeper with | cos 8* | for resolved processes than for direct ones. 



As discussed in the introduction and in section |3.1.3| , with a cut in x° bs it is possible to 



make a separation between events that are likely to be direct from those that are likely to 
be resolved events. An example of this is shown in Fig. ^ where direct and resolved events 
generated by Pythia is compared with results presented by the ZEUS Collaboration f31 



with a good agreement in the shape of the distributions. The plots were produced by using 
the HzTool package. 

In the analysis of the ZEUS Collaboration a cut x° bs < 0.75 was used to classify the 
photon as resolved while for x° bs > 0.75 it was classified as a direct photon. Of the 
events classified as having a direct photon a significant contribution comes from events 
with an anomalous fluctuation of the photon. This is well understood since, in this 
component, the quark (antiquark) often carries a large fraction of the original photon 
momentum. The two jets are then likely to be produced at small scattering angles 6* 
as described by the matrix elements with gluon exchange. Using the SaS ID parton 
distribution, the anomalous component is important for either classification, whereas the 
VMD part is highly suppressed and negligible by the x° bs > 0.75 cut. For SaS 2D, the 
VMD contribution is small for events classified as direct but of equal importance as the 
anomalous part for x° bs < 0.75. The choice of scale \ii is not crucial at this low Q 2 values, 
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Figure 9: The dijet angular cross section dcr/d| cos#*| for resolved (x° bs < 0.75) and 
direct (x° bs > 0.75) processes. 9* is the jet scattering angle w.r.t. the beam axis 
measured in the dijet centre of mass system. Q 2 < 4 GeV 2 and 0.25 < y < 0.8. 
The SaS ID and SaS 2D curves give the total contribution from the resolved and 
direct components; the VMD and direct components are also shown separately. The 
anomalous contribution can be obtained by subtracting the direct and the VMD 
contribution from the full one. 

the /x 3 scale is used and the results with /x 5 are within errors. 

The possibility for the incoming electron to emit extra photons, that may go unde- 
tected, is not taken into account by the model. This effect decreases the lepton energy; 
hence the photon energy calculated here would be overestimated. A simple estimation of 
these effects give corrections to the photon energy at the percent level. The consequence 
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would be that the given x° bs distribution will be shifted towards higher values. An in- 
clusion of such effects would shuffle some of our resolved events into the direct category, 
thus improving the description of the relative amounts. This could be solved by using 
lepton-inside-lepton structure functions, though the implementation of it is postponed for 
the future. 

From these observations it is tempting to conclude that a cut in x° bs does not so 
much separate events with the resolved anomalous component (perturbatively calculable 
fluctuation) of the photon from the direct component, but rather separates it from the 
VMD component (non-perturbative fluctuation) of the photon. This statement is true to 
leading order but the anomalous fluctuation contain pieces that a next-to-leading-order 
calculation of the direct component would give. Hence, in this sense also anomalous events 
can be considered as a contribution from the direct component making the x° bs well suited 
to separate the direct from the resolved component. Nevertheless, the classification of the 
photon into different components used throughout in this paper is based on the leading 
order description, viewing the anomalous fluctuation as a resolved photon. 



3.2.2 Inclusive ep Jet Cross Sections 

Inclusive ep jet cross sections have been measured by the HI collaboration in the 
kinematical range < Q 2 < 49 GeV 2 and 0.3 < y < 0.6. The differential jet cross 



sections da ep /dE]_ } da ep /dri* and the inclusive 7*p jet cross section in Fig. [K], |TT| and |12 
respectively, were produced with the Hz Tool package. The E]_ and rj* are calculated in 
the 7*p centre of mass frame where the incident proton direction corresponds to positive 
r]*. 

For dcr e p/d-E^L and d<7 ep /d?7* data is available in nine different Q 2 bins, four of them are 
shown here with similar results for the intermediate bins. The SaS ID parton distribution 
together with a few different Hi scales; Hii /i 3 , /i 5 and /z 6 , are used to model the resolved 
photon component. The other choices of scales, H2 and /14, interpolates between these 
results. 

In the highest Q 2 bin the direct component is the dominant contribution; the virtuality 
of the photon is for most events of the order of or larger than the transverse momenta 
squared, Q 2 However, the resolved component is not negligible and all the scales 

Hi, except Hi, depend on the photon virtuality. This gives a larger resolved component in 
this region as compared to the the conventional choice, Hi = P±- m t ne low-Q 2 bin the 
Hi scales do not differ much from p±, i.e. the results are not sensitive to the scale choice 
Hi- The exception is He, which there overshoots the data. The h\ = P± + Q 2 /^ scale (not 
shown) gives nice agreement with data for all different Q 2 bins. 

Since the VMD part dies out quickly with increasing photon virtuality, multiple inter- 
actions will only be visible at low Q 2 (multiple interactions for the anomalous component 
is not in the model so far). The anomalous component dominates over the VMD com- 
ponent already at 1 GeV 2 . Therefore, multiple interactions for the VMD component can 
safely be neglected for the distributions shown in this section. 

Changing the photon parton distribution from SaS ID to SaS 2D will give a slightly 
lower result for the low-Q 2 bins, Fig. O. The GRS LO |22| parton distribution gives a 



similar result as the other two. It has a restricted regime of validity, Q <C h (imple- 
mented as 5Q 2 < h 2 ) and Q 2 < 10 GeV 2 . Therefore, it is absent in the high-Q 2 bin, and 
some regions of phase space are cut out in the low--E^ bins and in the rj* distributions 
(5Q 2 ^ (E*[) 2 ). One could imagine larger differences for virtual-photon parton distribu- 
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Figure 10: The differential jet cross section da ep /dE^_ for jets with —2.5 < rj* < —0.5 
and 0.3 < y < 0.6. 



tions that from the onset are more different, so the good SaS/GRS agreement is somewhat 
fortuitous. Using a parton distribution for a real photon cannot describe the Q 2 damp- 
ening in the distributions shown in this section. As an example of this the GRV LO [Q 
distribution has been used; the change here is solely from the photon flux. 

Using CTEQ 3L instead of GRV 94 LO as the proton parton distribution reduces the 
result in some E* ± and rf bins by half, Fig. [TJ]. The GRV 94 HO parton distribution give 
a slightly lower result (as compared to GRV 94 LO). The differences mainly come from 
the gluon distributions, that are not yet so well constrained from data. In the modeling 
of the parton distributions, it is a deceptive accident that the more well-known proton 
parton distribution gives a larger uncertainty than the photon one. It offers a simple 
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Figure 11: The differential jet cross section da ep /drj* for jets with E]_ > 5 GeV and 
0.3 < y < 0.6. 



example that also phenomenology of other areas may directly influence the interpretation 
of photon data. 

3.2.3 Inclusive 7*7* Jet Cross Sections 

The OPAL collaboration has measured inclusive one-jet and two-jet cross sections in the 
range |?^ e1: | < 1 and requiring E 3 ® to be larger than 3 GeV [35]. The centre of mass 



energies were 130 and 136 GeV. In the analysis, a cone jet algorithm was used with a 
cone radius R = 1 and -Ej 6 ^ = 2 GeV; differences in the application of this algorithm 
may affect the results (we used the PXCONE algorithm from HzTool). The two-jet cross 
sections were obtained by measuring events with at least two jets, and then using only 
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Figure 12: The inclusive 7*p jet cross section cr 7 . p for jets with —2.5 < rj* < —0.5 
and 0.3 < y < 0.6. 

the two hardest jets. In Figs. |I5|-|T8|, the inclusive jet cross sections as a function of E 3 ± 



or rf et are shown, with events generated at y/s^ = 133 GeV. The anti-tagging conditions 
at this energy corresponds to a Q^ax = 0-8 GeV 2 . Due to migration effects p^ r ^ was 
chosen to 1.5 GeV, and p± l min = 1.4 GeV. The symmetric cuts used gives no difference 
between the two single-resolved contributions, instead the direct- VMD (VMD-direct) 
and direct-anomalous (anomalous-direct) contributions is shown together with the direct 
and double-resolved ones. 

At low E 3 ^ the double-resolved events are dominating and at larger E 3 ± it is the 
direct processes since more energy goes into the hard scattering in the latter case. For 
single-resolved events, the SaS ID VMD component dies out much quicker with increas- 
ing; E?f than the SaS 2D one which is comparable with the direct-anomalous events at 
high E 3 ^". For both cases, at high E 3 ^, the direct-anomalous components give the same 
order of magnitude contribution to the cross section as the double-resolved events. The 
biggest difference between the two parton distributions can be seen at low E 3 * and for 
the \rf et \ distributions, where the double-resolved events dominates; it is a reflection of 
the difference in normalization among the contributions. For the SaS 2D case, this kine- 
matical region makes the VMD component more important than the anomalous one; as 
a consequence multiple interactions play an important role. The double-resolved contri- 
bution for SaS 2D without multiple interaction is half of the one shown here. Clearly, 
for the SaS ID case the opposite is true: the importance of the components are reversed. 
In the region of high E 3 ^, where the direct events dominate, the model is undershooting 
data. On the other hand, there is nice agreement with data for the \rf et \ distribution when 
using SaS 2D. 
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Figure 13: The differential jet cross sections da ep /dE]_ and da ep /drj* for jets with 
—2.5 < r]* < —0.5 and 0.3 < y < 0.6. Only resolved events are used to show the 
sensitivity to the photon parton distribution. 



3.2.4 Forward Jets in ep 

Jet cross sections as a function of Bjorken-x, x-Qj, for forward jet production (in the proton 
direction) have been measured at HERA |36|. The objective is to probe the dynamics 



of the QCD cascade at small XBy The forward jet is restricted in polar angle w.r.t. the 
proton and the transverse momenta should be of the same order as the virtuality of 
the photon, suppressing an evolution in transverse momenta. If the jet has a large energy 
fraction of the proton, there will be a big difference in x between the jet and the photon 
vertex, x-q] <C Xj Ct , allowing an evolution in x. The above restrictions will not eliminate 
the possibility of having a resolved photon, although the large Q 2 values are not in favour 



25 



0.65 <Q 2 < 1.2 GeV 2 



36 < Q 2 < 49 GeV 2 



le+00 



> 
O 

■g le-01 

* H 



"D 



HI — •- 

GRV94LO 

GRV94HO 

CTEQ3L 



le-01 



le-02 



le-02 



le-03 



le-03 



4 6 8 10 12 14 16 18 20 
Ex (GeV) 




10 12 14 16 18 20 
Et (GeV) 



le+00 



"D 



le-01 



le-01 



hi f 



-2.5 



-1.5 

* 

11 



-0.5 



-2.5 



-1.5 

* 



-0.5 



Figure 14: The differential jet cross sections da ep /dE* ± and da ep /drj* for jets with 
—2.5 < rj* < —0.5 and 0.3 < y < 0.6. Only direct events are used to show the 
sensitivity to the proton parton distribution. 



of it. 

The HzTool routines were used to obtain the results in Fig. 19. Five different scales 
are shown. A larger forward jet cross section is obtained with a stronger Q 2 dependence, 
with the scale //§ = p\ + Q 2 in best agreement with data, in general agreement with the 
conclusions in |24|]0. The constraint (p 3 ±) 2 — Q 2 gives a large enhancement for the \i\ 
and //5 scales, which are combinations of (pP arton ) 2 an d Q 2 , whereas \x<i and \i$ have a Q 2 

3 However, we note a not complete agreement in the description of the inclusive jet cross sections. 
In |p4| a different option of the SaS extension to virtual photons was used than here, but we have checked 
that this cannot be the main reason for the discrepancy. 
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Figure 16: The inclusive two-jet cross section as a function of E 3 ^_ ; for more infor- 
mation see Fig. |T5| . 

dependence that are scaled down by s. The choice of scale does not only affect the resolved 
photon contribution but also the direct photon, arising from the scale dependence in the 
proton parton distribution, as seen in Fig. ^0[ As a check for the direct component, a 
simple comparison was made between Pythia, using the /zi = p± scale and the direct 
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Figure 18: The inclusive two-jet cross section as a function of rf et ; for more infor- 



mation see Fig. 17. 



processes, and LEPTO 0, without soft colour interactions, giving similar results. The 
rather large Q 2 values, Q 2 ~ (p^ ) 2 , suppresses VMD photons and favours the SaS ID 
distribution which is the one used here, though the difference is small. 

Note that the scale undershoots the forward jet cross section data and overshoots 
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Figure 19: Forward jet cross section as a function of x compared with HI data (with 
statistical and systematic uncertainties added in quadrature). Five different scales 
are shown at two different p 1 ^ cuts, 3.5 and 5 GeV. x- iCt > 0.035, 0.5 < (p^ ) 2 /Q 2 < 2 



and 7° < jet < 20°. 
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Figure 20: Same as in Fig. 19 but only the direct component is compared with data 



for two different scale choices. 

the inclusive jet distributions at low Q 2 , so it is not a real alternative. As a further 
check, with more data accumulated and analyzed, the (j?±) 2 /Q 2 interval could be split 
into several subranges which hopefully would help to discriminate between scale choices. 
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3.2.5 Importance of Longitudinal Resolved Photons 

In this section we will study the importance of longitudinal resolved photons. The differ- 
ent .R-factors described in section |2]3| will be used to provide some estimates. Of those 
distributions studied so far, we will concentrate on the da cp /dE*[ distributions for inclu- 
sive jets and the d<T ep /dx for forward jet cross sections. A sensible Q 2 -dependent scale 
choice, [A3, together with the SaS ID distribution will be used throughout. 

With a = 1 the different alternatives are shown in Fig. |2l] for the dcr ep /d^ distribu- 
tions together with the result from pure transverse photons, i.e. a = 0. The importance 
of the resolved contributions decreases with increasing Q 2 , see Fig. [10], which makes the 
asymptotic behaviour less crucial. The onset of longitudinal photons governed by the R\ 
and i?2 alternatives are favoured whereas the R3 one overshoots data in the context of 
the other model choices made here. 

In Fig. ^ the same alternatives are shown for the forward jet cross sections. With 
this scale choice, /i 3 , none of the longitudinal resolved components (together with the 
direct contribution) are sufficient to describe the forward jet cross section. The resolved 
contribution with R 3 is about the same as the one obtained with the scale /x| = p\ + Q 2 
(without longitudinal contribution); the difference in the total results originates from the 



difference in the direct contributions, see Fig. ^0[ With Ri and a = 1, the /i 5 scale (not 
shown) overshoots the data, but undershoots in combination with /x 4 . The \x± scale in 
combination with i? 3 and a = 1 is in nice agreement with data. 

The above study indicates, as expected, that longitudinal resolved photons are im- 
portant for detailed descriptions of various distributions. It cannot by itself explain the 
forward jet cross section, but may give a significant contribution. Combined with other 
effects, for example, different scale choices, parton distributions, underlying events, etc. 
it could give a reasonable description. The model (s) so far does not take into account 
the difference in x distribution or the k 2 scale (of the 7* — > qq fluctuations) between 
transverse and longitudinal photons. As long as the distributions under study allow a 
large interval in x the average description may be reasonable. In a more sophisticated 
treatment these aspects have to be considered in more detail. 

3.3 Example of Further Tests: Forward Jets in e + e 

With the experience of forward jets at HERA, we suggest a similar study at LEP. The 
optimal kinematical and forward jet constraints have to be set by each collaboration 
itself; the study here will give the order of magnitude for the cross section and point out 
uncertainties in the model. 

Comparing with forward jets at HERA, one of the leptons will play the role of the 
proton. Some of the constraints will be taken over directly, kinematic cuts: 

y > 0.1, E' e > 11 GeV, (46) 



and jet selection: 



> 3.5 GeV (47) 

> 0.035 (48) 

0.5 < ^-<2 (49) 

Q 

3° < 9 ]ct < 20° (50) 
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Figure 21: The differential jet cross section da ep /dE^_ for jets with —2.5 < rj* < —0.5 
and 0.3 < y < 0.6. 



To fulfill the jet selection one of the leptons has to be tagged in order to know the virtuality 
of the photon. To obtain a reasonable number of events the other lepton is not tagged, 
imposing a here chosen to 1.5 GeV 2 . With a centre of mass energy of 200 GeV, the 

smallest accessible xb] — ^ is around 10~ 4 , where Q 2 and y is calculated from the tagged 
electron, omitting the virtuality of the other photon. In a more sophisticated treatment 
also double-tagged events are analyzed; then one of the photons plays the role of a proton 
and the forward jet should be defined with respect to one of the photons. 

A cone jet algorithm with cone radius R = 1 and ffif m - m = 2 GeV is used for jet finding. 
As for the case at HERA, the /i 5 scale gives the largest forward jet cross section, about 
twice as large as with the /i2 scale, Fig. |2^. Most of the differences arise from the double- 
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Figure 22: Forward jet cross section as a function of x. 



resolved events. Double-resolved and single-resolved events, where the resolved photon 
give rise to the forward jet, dominate the forward jet cross section, Fig. E3| and |24l. At 



low x, for the n§ scale, the double-resolved contribution is close to an order of magnitude 
larger than the direct one. For the fi2 scale it is about a factor of four. As for the case 
at HERA, the rise of the forward jet cross section at small x is dominated by resolved 
photons. A study like this at LEP could be an important cross check for the understanding 
of resolved photons and that of small- x dynamics. 



4 Summary and Outlook 

The field of photon physics is rapidly expanding, not least by the impact of new data 
from HERA and LEP. The prospects of building a Linear Collider, with its objective of 
high-precision measurements and to search for possible new physics, requires an accurate 
description of photon processes. The plan here is to have a complete description of the 
main physics aspects in 7p and 77 collisions, which will allow important cross checks to 
test universality of certain model assumptions. As a step forward, we have in this article 
concentrated on those that are of importance for the production of jets by virtual photons, 
and are absent in the real-photon case. While we believe in the basic machinery developed 
and presented here, we have to acknowledge the many unknowns — scale choices, parton 
distribution sets (also those of the proton), longitudinal contributions, underlying events, 
etc. — that all give non-negligible effects. To make a simultaneous detailed tuning of all 
these aspects was not the aim here, but rather to point out model dependences that arise 
from a virtual photon. 

When Q 2 is not small, naively only the direct component needs to be treated, but in 
practice a rather large contribution arises from resolved photons. For example, for high-Q 2 
studies like forward jet cross sections, Fig. p^pO|, or inclusive differential jet cross sections, 
Fig. |T0|-|r2]. Resolved longitudinal photons are poorly understood and the model presented 
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Figure 23: Forward jet cross section as a function of x. p*± > 3.5 GeV, Xj et > 0.035, 
0.5 < )7Q 2 < 2 and 3° < 6> jet < 20°. 

here can be used to estimate their importance and get a reasonable global description. 
Longitudinal effects are in most cases small but of importance for fine-tuning. 

In the study of the dijet angular distributions in photoproduction the relative amount 
of direct and resolved events is not so well described by the model. In the future, this 
could be improved by using lepton-inside-lepton structure functions. 

The inclusive 7*7* one-jet and two-jet cross sections are well described except for the 
high E 3 ± region of the E 3 ^ distribution. In this region, the direct events are dominating. 
Currently, owing to the lesser flexibility in the modeling of the direct component, we do 
not see any simple way to improve the model. The factorized ansatz made for the photon 
flux is expected to be valid in this kinematical range; interference terms are suppressed 
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Figure 24: Forward jet cross section as a function of x. tff" > 5 GeV, .r 
0.5 < {^f/Q 2 < 2 and 3° < jet < 20°. 



jet > 0.035, 



by Q\Q\jW 2 *^ • However, differences in the application of the cone jet algorithm may 
affect the results. 

The forward jet cross section presented by HI |5(J is well described by an ordinary 
parton shower prescription including the possibility of having resolved photons. The 
criteria that the p 3 ^ should be of the same order as Q 2 , makes the scale choice crucial and 
A 4 ! = v\ + Q 2 is favoured by data, as concluded in []24] . With this experience we predict 
the forward jet cross section to be obtained at LEP. With more data accumulated and 
analyzed, the (p 3 ±) 2 /Q 2 interval could be split into several subranges, which hopefully 
would help to discriminate between different scale choices. 

Multiple interactions for the anomalous component are not yet included, and is not 
expected to be of same importance as in the VMD case. However, for low k 2 fluctuations 
it may be important, especially for SaS ID, and need to be investigated. 

After this study of jet production by virtual photons it is natural to extend the model- 
ing to \ow-pj_ events. Clearly, a smooth transition from perturbative to non-perturbative 
physics is wanted. One idea is to make use of a parameterization of the total 7p and 77 
cross section in terms of a pomeron and a reggeon exchange. Starting from the real-photon 
case, dipole dampening factors are introduced for the generalization to virtual photons. 
For example, the 7p cross section is divided into a VMD, an anomalous, a direct and a 
DIS 7*q — > q process part. In the limit Q 2 —>■ 0, the first three event classes remain. On 
the contrary, when Q 2 increases from zero to high Q 2 ; the resolved processes dies out (as 
given by the dipole factors), the direct also drops and finally only the DIS process remain. 
At intermediate Q 2 values, the direct processes and the DIS (+parton showers) process 
overlap, since, in some regions of phase space, they are equally valid descriptions of the 
same physics. It thus becomes necessary to avoid double-counting, e.g. by introducing 
Sudakov style form factors for the DIS process, suppressing those parton configurations 
covered by the direct processes. We intend to return to this issue in a future publication. 
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